arXiv: 1507.02565v2 [hep-ph] 2 Sep 2015 


ZU-TH 20/15 
MITP/15-051 


Transverse-momentum resummation 
for vector-boson pair production at NNLL-j-NNLO 


Massimiliano Grazzini^Q Stefan Kallweit^, 
Dirk Rathlev^) and Marius Wiesemann^ 


( a ) Physik-Institut, Universitat Zurich, CH-8057 Zurich, Switzerland 

^ PRISMA Cluster of Excellence, Institute of Physics, 
Johannes Gutenberg University, D-55099 Mainz, Germany 


Abstract 

We consider the transverse-momentum (px) distribution of ZZ and W + W~ boson 
pairs produced in hadron collisions. At small px, the logarithmically enhanced 
contributions due to multiple soft-gluon emission are resummed to all orders in QCD 
perturbation theory. At intermediate and large values of px , we consistently combine 
resummation with the known fixed-order results. We exploit the most advanced 
perturbative information that is available at present: next-to-next-to-leading logarithmic 
resummation combined with the next-to-next-to-leading fixed-order calculation. After 
integration over px, we recover the known next-to-next-to-leading order result for the 
inclusive cross section. We present numerical results at the LHC, together with an 
estimate of the corresponding uncertainties. We also study the rapidity dependence of 
the px spectrum and we consider px efficiencies at different orders of resummed and 
fixed-order perturbation theory. 
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1 Introduction 


Run 1 of the Large Hadron Collider (LHC) has been a great success for the Standard Model (SM). 
The collected data are in good agreement with the theoretical predictions so far and led to the 
discovery |1,2| of a resonance at a mass of 125 GeV, which appears to be fully consistent with 
the SM Higgs boson. Among the most important reactions at hadron colliders is the production 
of vector-boson pairs. This class of processes gives access to the vector-boson trilinear couplings 
which may be modified in a large set of Beyond the Standard Model (BSM) theories. Even small 
deviations in both the production rate and the shape of distributions could be a signal of new 
physics. Anomalous couplings related to vector-boson pair production have been constrained first 
by LEP2, and later by the Tevatron for larger invariant masses. ATLAS and CMS will continue to 
tighten the bounds on anomalous couplings, especially with increasing sensitivity during Run 2 of 
the LHC0 


On the other hand, vector-boson pair production constitutes an irreducible background to 
new-physics searches as well as Higgs studies. Particularly important are the off-shell effects below 
the ZZ and W + W~ thresholds, relevant for the Higgs signal region, and the high-mass tail used to 
extract the width of the Higgs boson |4[j6]. Furthermore, Higgs boson measurements, in particular 
in the W + W~ channel, strongly rely on the background rejection through specific categories based 
on the transverse momenta of final-state particles, such as the classification into jet bins or in 
the Higgs transverse momentum and related variables. An accurate modelling of the respective 
observables for both signal and backgrounds is crucial for such analyses. 


The first precise predictions for ZZ production at hadron colliders in the SM were obtained at 
the next-to-leading order (NLO) already more than 20 years ago for stable Z bosons 7,8j, and the 
leptonic decays were added in Ref. j9j. The full spin correlations and off-shell effects at the NLO 
were first included in Refs. 10,11 using the corresponding one-loop helicity amplitudes Il2|. An 


important loop-induced contribution proceeds through gluon fusion; it is enhanced by the gluon 
densities and was first computed for on-shell Z bosons in Refs. 13,14 , while leptonic decays were 
included later 15 1T|. All the contributions to ZZ production discussed so far are implemented in 
the numerical program MCFM 18 . Electroweak (EW) corrections were evaluated in Ref. 19,20 


while on-shcll ZZ+ jet production is known through NLO QCD 21,22 . Recently, the inclusive 


cross section for the production of ZZ at the next-to-next-to-leading order (NNLO) has been 
presented 


23 


The production of W + W~ pairs constitutes the largest cross section among all massive vector- 
boson pair production modes. On the other hand, its leptonic decay (W + W~ —>- l + l~vv) embodies 
the most challenging experimental signature since the presence of two neutrinos prohibits a 
reconstruction of mass peaks. Therefore, a precise understanding of both the signal and the 
background is required. 


The W+W 
component 


cross sections for on-shell W bosons at the NLO 24 25 as well as the gluon-fusion 


13 


have been known for decades. Also in this case, spin and off-shell effects were 
included in the NLO prediction 10,11 after the relevant one-loop helicity amplitudes had been 
computed 12 . Leptonic decays for the loop-induced gluon-fusion contribution were considered 


in Refs. More recently, also interference effects with the Higgs production mode through 

tSee Ref. 0 and references therein. 
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gluon fusion were determined 29 


production discussed so far are implemented in MCFM 

On-shell W+W 


Analogously to ZZ production, all the contributions to W+W 

Furthermore, EW corrections have 


been evaluated 20 30,31 


18 


through NLO QCD in Refs. 32-34 


association with up to one jet at NLO have been presented in Ref. 
results for the inclusive W+W~ cross section have been obtained 


production in association with one jet has been studied 
Detailed Monte Carlo simulations of e+v e prv^ production in 

Recently, the first NNLO 


35 


36 


Due to the very recent computation of the qq —» VV' helicity amplitudes at two-loop order 
the inclusion of the off-shell effects and the leptonic decays in the NNLO cross section is 


37,38 


expected in the near future. In the meanwhile, also the calculation of gg —> VV' helicity amplitudes 
has been performed at two-loop order 39 401. This renders the evaluation of NLO QCD corrections 


to the gluon-fusion channel feasible. 

The production of ZZ pairs in hadron collisions has been measured extensively at the Tevatron 
and the LHC (see Refs. 41-47 for some recent results). The W+W~ cross section has also been 
measured already at the Tevatron, see e.g. Ref. |48], and at the LHC both at 7 TeV |49,50 and 8 
TeV [46,|5llj52|. The ATLAS collaboration recently reported an excess 51 with respect to the SM 


prediction, which has drawn a lot of attention on the W+W process, since the W+W final state 

In the meanwhile, the excess has been alleviated 


is a typical signature in many BSM scenarios 53 


to a significant degree by the recent NNLO computation 36 
the CMS collaboration 


The more recent measurement by 


52 is in good agreement with the NNLO prediction. 


The transverse-momentum distributions of ZZ and W+W pairs are among the most important 
differential observables for these processes. The px spectrum has already been measured in the 


case of ZZ production 47 at the LHC. Transverse-momentum resummation for ZZ and W+W 
production has been studied in Refs. 54-58 . In all these calculations the resummed computation 


is essentially performed up to next-to-leading logarithmic (NLL) accuracy (next-to-next-to-leading 
logarithmic (NNLL) effects in the Sudakov exponent are considered in Refs. 155, 57,;58]) and 
matched to the fixed-order result up to the first order in the QCD coupling og. 

In this paper we consider transverse-momentum resummation for the production of both ZZ 
and W+W~ pairs. We use the formalism of Ref. [591 to perform the first computation of the 
Pt spectrum at full NNLL accuracy matched to the d(a|) fixed-order result valid at large px- 
Although our focus is on the inclusive px spectrum of the ZZ and W+W~ system, our computation 
is fully differential in the degrees of freedom of the vector bosons and allows us to include their 
EW decays, once the helicity amplitudes are implemented The r es u mrri at ion formalism of 
Ref. (59] is closely related to the subtraction method of Ref. |63|, which was used to compute the 
NNLO cross section for these processes 23,361. For W+W~ production we employ the four-flavour 


scheme (4FS) and remove all contributions with final-state bottom quarks from our computation 
of the W+W~ transverse-momentum distribution in order to eliminate the contamination from 
tt and Wt production. The difference to the prediction in the five-flavour scheme (5FS), where 
such terms have to be consistently subtracted, has been shown to be small for the NNLO inclusive 
cross section 36 . Furthermore, we neglect the loop-induced gluon-fusion contribution throughout 


this paper, since, up to NNLO, it contributes only at px = 0. 


We note that in the CMS measurement reported in Ref. 52 an approximate NNLL prediction 


Mhe analogous computations for Higgs, single vector-boson production, and diphoton production are presented 


in Refs. 60 , 61 and 62 , respectively. 
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58 of the pt spectrum has been used to correct the spectrum from the Monte Carlo simulation. 


Along these lines, the computation reported in this paper will be useful in order to validate the 
predictions obtained from Monte Carlo simulations for both ZZ and W + W~ production, as done 
in the case of Higgs boson production with the calculation of Ref. 


59 


The manuscript is organized as follows. In Section [2] we review the transverse-momentum 
resummation formalism applied to vector-boson pair production. In Section [3] we report our 

In Section 13.21 we present our numerical predictions for the inclusive pr spectrum and study 
the ensuing uncertainties. In Section [3731 we analyse the behaviour of the spectrum at different 
rapidities of the vector-boson pair. In Section l3~4l we investigate pt efficiencies at different orders 
in resummed and fixed-order perturbation theory. In Section |4] we summarize our results. 


numerical results, starting with remarks on the choice of the resummation scale in Section 3.1 


2 Transverse-momentum resummation for vector-boson 
pair production 


In this Section we recall the main points of the transverse-momentum resummation formalism we 
use in this paper. For a more detailed discussion the reader is referred to Refs. 59,64,65 . 


We consider the inclusive hard-scattering process 


hi(Pi) + h 2 (P 2 ) V(p 3 ) + V\p A ) + X , 


( 1 ) 


where the collision of the two hadrons h\ and h 2 with momenta P\ and P 2 produces the two vector 
bosons of momenta p 3 and p±. In the center-of-mass frame the momentum of the vector-boson pair 
q = p 3 + P4 is fully specified by the invariant mass M 2 = (p 3 + p 4 ) 2 , the rapidity y = \ In ^||, and 
the transverse-momentum vector p T . 


The kinematics of the vector bosons is fully determined by the vector-boson pair momentum 
g M = p 3 +P 4 (with p\ = rriy and p\ = rriy,), and by the additional and independent variables that 
specify the angular distribution of the vector bosons with respect to Throughout this paper we 
always consider quantities which are inclusive over these angular variables. 


According to the QCD factorization theorem, the differential cross section can be written as 


da vv ' f 1 f 1 

dM 2 dp 2 dy (S/» Pr,M,s) =Y^ l o dx l J o dx 2 f ai/hl (xi, V 2 f) fa 2 /h 2 (x2, V 2 f) 


X 




dM 2 dpy dy 


(y,p T ,M,s,a s (p 2 R ),p%p 2 F ), 


( 2 ) 


where f a /h(x, Pf) (a = q, q, g) are the density functions of parton a in hadron h at the factorization 
scale hf'i Pr is the renormalization scal^j da^' 2 is the partonic cross section. The rapidity y and 
the center-of-mass energy s of the partonic scattering process are related to the corresponding 

§ Throughout the paper we use parton densities in the MS factorization scheme and as{q 2 ) is the QCD running 
coupling in the MS renormalization scheme. 


3 
















hadronic variables y and s by 


1 , Xi 

y = y- o ln ~> 

2 x 2 


S = X\X 2 s , 


(3) 


When the transverse momentum px of the vector-boson pair is of the same order as the invariant 
mass M, the QCD perturbative expansion is controlled by a single expansion parameter, «g(M), 
and fixed-order calculations can be safely applied. In this region, QCD radiative corrections are 
known to O(o§) 121,22,32-34 . When px M the convergence of the perturbative expansion is 


spoiled by the presence of large logarithmic terms of the form aig ln m (M 2 /p^), that need to be 
resummed to all orders. 


The resummation is performed at the level of the partonic cross section, which is decomposed 


as 


d*Z> 


j ~VV' ,(res.) ,.W',(fin.) 

LLU a±a2 ^ UU aiCL2 


_ ('ll 

dM 2 dp\dy dM 2 dpxdy dM 2 dp\dy 

The first term on the right-hand side of Eq. Q contains all the logarithmically-enhanced contribu¬ 
tions at small px and has to be evaluated by resumming them to all orders. The second term is 
instead free of such contributions and can thus be evaluated at fixed order in perturbation theory. 

Resummation is based on the factorization of soft and collinear radiation and is viable in the 
impact-parameter ( 6 ) space, where the kinematical constraint of momentum conservation and the 


factorization of the phase space can be consistently taken into account 66-68 . Using the Bessel 


transformation between the conjugate variables px and b, the resummed component is expressed 
as 


59,68 


ipVV',(ves.) 114-2 roo h 

dmtpIdy = ~T l db 2 MbpT) S. M, i; as, A, A) . 


(5) 


where Jo(x) is the 0-order Bessel function. In the case of fully inclusive px resummation, the rapidity 
dependence is integrated out in Eq. ([ 5 ]). I 11 that case it is convenient to consider Mellin moments 
with respect to the variable z = M 2 /s. However, in order to retain the rapidity dependence in the 
resummed cross section we take the ‘double’ (N\, N 2 ) Mellin moments with respect to the variables 
z\ = e + ^M/\/§ and z 2 = e~^M/\/~§ at fixed M, 


Wj^ 2 ) ( 6 ,M;a s , 4 , 4 )= / dz lZ ^~ l dz 2 z^~ 1 W 


w 


(b, y, M, s; as, dw dV ), ( 6 ) 


and organize the structure of Wyy/ in the following exponential form [64], 

a S: d 2 Ri /^f) = ^L{Ni,N2) {M'j a S; M 2 / : p 2 j, M 2 /pp, M 2 / Q 2 ) 
x exp{£ (JVl ,jv 2 )(a s , L, M 2 / p 2 R , M 2 /Q 2 )} , 

where we have defined the logarithmic expansion parameter L as 

Q 2 b 2 


L — In 


b 2 

u o 


(7) 


( 8 ) 


and b 0 = 2e lE ( 7 e = 0.5772... is the Euler number). The scale Q appearing in Eqs. ®§ , named 
resummation scale in Ref. 59 , parametrizes the arbitrariness in the resummation procedure, and 


has to be chosen of the order of the hard scale M. Variations of Q around its reference value 
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can be exploited to estimate the size of yet uncalculated higher-order logarithmic contributions. 
Therefore, Q plays a role very similar to fip and /xr for missing perturbative terms. The function 
does not depend on the impact parameter 6, and therefore includes all the perturbative 
terms that behave as constants as b —* oo. It can thus be expanded in powers of as = as(g R )- 


n 


■{N u N 2 )\ 


+ 




( 9 ) 


M z /Q z 


where a VV '^A is the partonic leading-order (LO) cross section. The exponent Q{n 1 ,n 2 ) includes the 
complete dependence on b and, in particular, it contains all the terms that order-by-order in as 
are logarithmically divergent as b — y oo. The logarithmic expansion of Gn reads 

G {Nl ,N 2 )(as(fi 2 R ), L- M 2 //4 M 2 /Q 2 ) = LgW(a s L) + g^ N 2 ) (a s L- M 2 /&, M 2 /Q 2 ) 

+ ~9 , (^i,Ar 2 )( Q! S-^’ -^ 2 /VL M 2 /Q 2 ) + . .. , (10) 


where the term L gA') collects the LL contributions, the function g^ includes the NLL contri- 
(3) 

butions, g^ controls the NNLL terms and so forth. 


The resummation of the large logarithmic terms carried out in Eq. ([7]), after transforming 
back to Pt space, allows us to obtain a well behaved transverse-momentum spectrum as px —> 0. 
However, the logarithmic expansion parameter L in Eq. (|8]) is divergent as b —» 0. This implies 
that the resummation produces higher-order contributions also in the high-p^ region, which is 
conjugated to b —y 0 after Fourier transformation. In this region the fixed-order cross section is 
perfectly viable and any resummation effect is necessarily artificial. To reduce the impact of such 
contributions, the logarithmic variable L is replaced by 


59 


L L, 


L = In 



( 11 ) 


The variables L and L are equivalent when Qb 1 but they have a very different behaviour as 
6—^0. When Qb 1 , L — > 0 and G(N!,n 2 ) 1 - Moreover, since the behaviour of the Sudakov 

form factor at b = 0 is related to the integral over px, the replacement in Eq. ( fTTj ) allows us to 
enforce a unitarity constraint such that the fixed-order prediction is recovered upon integration 
over px ■ 


A well known property of the formalism of Ref. 59 is that the process dependence (as well as 
the factorization scale and scheme dependence) is fully encoded in the hard function %vv'■ In 
other words, the functions g^ are universal: they depend only on the channel in which the process 
occurs at Born level (qq annihilation in the case of vector-boson pair production). Their explicit 
expressions up to i = 3 are given in Ref. 59 in terms of the universal perturbative coefficients A^q\ 
Aq 2 \ Aq 3 \ In particular, the LL function r/ 1 -’ depends on the coefficient the NLL 


function g^ also depends on A^’ and Bq L> 69 and the NNLL function g [ (RT] also depends 


| ( 2 ) 


i(i) 


on Ag 3 - 1 


70 


and B 


( 2 ) 

q,N 


71 


73 


A3) 
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The hard coefficients TL 


VV' 


coefficients O0 ' ,(1 ' ) are known since long time 72, 73 


depend on the process we want to consider. The first order 

they can be obtained from the one-loop 
scattering amplitudes qq —> VV' by using a process independent relation. By exploiting the 

production, in Ref. 


and vector-boson 75 


65 


such a 


expressions of H for single Higgs 74 
relation has been extended to 0[a g): this implies that the two-loop amplitude for qq —> VV' is 
the only process-dependent information needed to obtain the coefficient dO v ’( 2 L 

We now turn to the finite component of the transverse-momentum spectrum, i.e. the second 
term on the right hand side of Eq. (4). Since do^ does not contain large logarithmic terms 
in the small-pr region, it can be evaluated by truncating the perturbative series at a given fixed 
order. In practice, the finite component is computed starting from the customary fixed-order 
(f.o.) perturbative truncation of the partonic cross section and subtracting the expansion of the 
resummed cross section in Eq. (J5]) at the same perturbative order: 


r 1 

LtO'0,10,2 


d*z: 


[ da VV’,(re s.) 1 

U'O' Oi<22 

dM 2 dp\dy 

f.o. 

dM 2 dp‘fdy 

f.o. 

dM 2 dp\dy 


( 12 ) 


J f.o. 


At least formally, this matching procedure between resummed and finite contributions guarantees to 
achieve a uniform theoretical accuracy over the entire range of transverse momenta. At large values 
of pt, the resummation procedure cannot improve the fixed-order result, and the resummation 
(and matching) procedure is eventually superseded by the customary fixed-order calculations. 

In summary, the inclusion of the functions g^\ N ^, 00 V 'V) i n the resummed component, 

together with the evaluation of the finite component at NLO (i.e. C?(as))j allows us to perform 
the resummation at NLL+NLO accuracy. This is the theoretical accuracy of the calculations of 
Refs. 55 


561. Including also the functions g^ and dO 1 '^ 2 \ together with the finite component 
i.e. 0(aD) leads to full NNLL+NNLO accuracy. Using the recently computed two-loop 


65 , we are now able to 


at NNLO 

amplitudes for qq —> VV', and the process independent relation of Ref. 
present the complete result for the transverse-momentum distribution of the vector-boson pair up 
to NNLL+NNLO accuracy. We point out that the NNLL+NNLO (NLL+NLO) result includes the 
full NNLO (NLO) perturbative contribution in the entire px range. In particular, the NNLO (NLO) 
result for the double differential cross section do/(dM 2 dy) is exactly recovered upon integration 
over px of the differential cross section at NNLL+NNLO (NLL+NLO) accuracy. 

We conclude this Section by adding a few comments on the way in which our calculation is 
actually performed. The practical implementation of Eq. (|4]) is done in the numerical program 
MATRIX^] which is an extension of the numerical program applied in the NNLO calculations of 
Refs. 


77 and based on a combination of the g-p-subtraction formalism 1631 with the 
MUNICti|jJcode. Since already performed within the g^-subtraction formalism, the extension of 
these calculations to compute the resummed cross section is conceptually quite straightforward, 
and is obtained by replacing the hard-collinear terms in the fixed-order computation by the proper 
all-order resummation formula of Eq. ([5]). This procedure is the same that was applied to perform 


the NLL+NLO calculations for W + W 55 and ZZ 56 production, and the NNLL+NNLO 


calculation for Higgs boson production of Ref. 1601. 

^Matrix is the abbreviation of “Munich Automates qT subtraction and Resummation to Integrate X-sections”, 
by M. Grazzini, S. Kallweit, D. Rathlev, M. Wiesemann. In preparation. 

II Munich is the abbreviation of “MUlti-chaNnel Integrator at Swiss (CH) precision”—an automated parton level 
NLO generator by S. Kallweit. In preparation. 
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To obtain the numerical results presented here, the resummed component of Eq. ([5]) is evaluated 
with an extension of the numerical program used for the calculation of Higgs production 1601, 
based on the earlier computations of Refs. 59,64 . The hard-collinear coefficients are obtained 


by exploiting the implementation of the corresponding virtual amplitudes for the production of 
on-shcll ZZ and W + W~ pairs in Ref. [23 and Ref. 36 , respectively, and the knowledge of the 


collinear coefficients relevant to quark-initiated processes [75 


The finite component of Eq. (12) is obtained from an NLO calculation of W'+jet, computed 


with the Munich code, which provides a fully automated implementation of the Catani-Seymour 
dipole formalism |80, 811 as well as an interface to the one-loop generator OpenLoops [821 
to obtain all required (spin and color-correlated) tree-level and one-loop amplitudes. For the 
numerically stable evaluation of tensor integrals we rely on the Collier library 83 , which is based 
on the Denner-Dittmaier reduction techniques 


84,85 


and the scalar integrals of [86] . To deal 
with problematic phase-space points, OpenLoops provides a rescue system using the quadruple- 
precision implementation of the OPP method in CutTools 87 , involving scalar integrals from 
OneLOop 


3 Results 


In this Section we present our results for the resummed transverse-momentum distributions of 
W + W~ and ZZ pairs. We compare our NNLL+NNLO predictions to the results at the NLL+NLO, 
and discuss the corresponding theoretical uncertainties. Additionally, we also study the rapidity 
dependence of the px cross section as well as the p^-veto efficiency. 

For the EW couplings we use the so-called scheme, where the input parameters are Gf, mw, 
rriz■ In particular we set Gf = 1.16639 x 10 -5 GeV' 2 , mw = 80.399 GeV, rriz = 91.1876 GeV. 
We use the NNPDF3.0 sets of parton distribution functions (PDFs) 1891 with oisirriz) = 0.118. 
At NLL+NLO and NNLL+NNLO the running of ag is evaluated at two- and three-loop order, 
respectively. For ZZ production we consider Nf = 5 massless quarks/antiquarks. For W + W~ 
production we make use of the 4FS, which allows us to split off all contributions related to 
bottom-quark final states in order to remove the tt and Wt contamination from our computation. 
This is straightforward in the 4FS, because the W + W~bb process is separately finite. 

We consider proton-proton collisions at yfs = 8 TeV. The central values of the factorization and 
renormalization scales are set to /ip = Hr = Po = 2 my. The choice of the central resummation 
scale Q 0 is discussed in the next subsection. 


3.1 Choice of the central resummation scale 


As discussed in Section (2j the resummation scale Q is the scale entering the large logarithmic 
terms we are resumming (see Eq. Q), and it plays the role of the scale up to which resummation 
is effective. In on-shell Higgs [59] and vector-boson 90 production, the scale is typically chosen 


*We note that an independent computation of these coefficients in the framework of Soft Collinear Effective 


Theory has been presented in Refs. 78 79 . 
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Figure 1: Invariant mass ( Mww ) distribution in W + W pair production at NLO (red, 
dashed) and NNLO (blue, solid). 


equal to half the mass of the heavy boson (i.e. Q = mu /2 for Higgs production and Q = my/ 2 in 
the case of single vector-boson production). Higher values of the scale lead to a worse matching 
at high pt■ The natural extension of this choice for vector-boson pair production is a dynamical 
resummation scale Q = Myy/2 , since Myy is the hardness of the process. This is indeed the 
choice that was adopted in the calculations of Refs. 


55,56 


The following considerations apply both to ZZ and W + W~ production, so we will focus on 
W + W~ production from now on. In Fig. [l] we consider the invariant mass distribution of the 
W + W~ pair at NLO and NNLO. We see that the distribution is strongly peaked in the threshold 
region, and that it quickly decreases as Mww increases. As a consequence, for most of the W + W~ 
events, M ww > 2mw- 


We can compare the transverse-momentum distributions obtained with a dynamical resum¬ 
mation scale Q = M W w /2, and a fixed resummation scale Q = mw- In Fig. [2] we show the ratio 
(blue, solid curve) of the Q = mw result over the Q = M ww / 2 result. The bands are obtained by 
varying the resummation scale around the central value by a factor of two. Considering the ratio 
of the central curves for px < 250 GeV, the differences between a fixed and a dynamical scale are 
extremely small and remain at the 1-2% level over the whole range. In this region of transverse 
momenta the uncertainty bands obtained with the two choices overlap and are similar in size. In 
fact, since Q = mw leads to slightly larger uncertainties, it appears to be the more conservative 
choice. Therefore, we can conclude that either choice of the resummation scale is perfectly valid 
and indeed consistent with each other as expected from the discussion of the invariant mass 
distribution. 


Looking further at the comparison of the high-pr tails in Fig.[2] (px > 250 GeV), we observe a 
very well known feature 59,60 90 92 of the applied matching procedure, namely the fact that 







































Figure 2: NNLL+NNLO transverse-momentum distribution of the W + W~ pair with 
a fixed scale Q = m\y normalized to the same cross section with a dynamical scale 
Q = M\yw /2. The bands are obtained by variation of the resummation scales in the 
numerator by a factor of two around the central scale. For reference, we show the 
fixed-order NNLO curve with the same normalization. 


for large values of the resummation scale the fixed-order cross section (black dotted curve) is 
not recovered in the tail of the distribution. It is important to recall that transverse-momentum 
resummation is supposed to improve the perturbative expansion in the low-pr region. At large pt, 
any large dependence on the resummation scale is necessarily artificial and an unwanted remnant 
of the matching procedure. This behaviour is precisely what we observe for the dynamical scale 
choice Q = Mww/% in Fig. [2] With this choice, in fact, the resummed result loses predictivity, as 
its uncertainty becomes increasingly large. By contrast, a fixed resummation scale Q = mw, which 
is always smaller than Q = M-^w/ 2, eventually leads to a more consistent high-pr behaviour of 
the resummed prediction. 

Based on the above results, we make Qq = my our default choice of the resummation scale in 
what follows. 


3.2 Inclusive transverse-momentum distribution 


We now present our resummed predictions for the inclusive transverse-momentum spectrum of the 
vector-boson pair and compare them with the corresponding fixed-order results. We concentrate on 
W + W~ production since we observe no saliently different features in the ZZ case. For completeness, 
we provide the corresponding reference prediction with uncertainties for ZZ below. 


Before presenting our resummed predictions, we recall the well known fixed-order results at 

In Fig.[3]we show the NLO and NNLO distributions together with their 


O(as) and 0(a g) |32 34 
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Figure 3: Transverse-momentum distribution of the W + W~ pair at NLO (red, dashed) 
and NNLO (blue, solid); thick lines: central scale choices; bands: scale uncertainty from 
Hf and Hr obtained as described in the text. Lower inset: results normalized to the 
NLO prediction at central values of the scales. 


perturbative uncertainties. The uncertainty bands are obtained by varying /if and p,R hr the range 
mw < {/if,/i,r} < 4 mw with the constraint 0.5 < Hf/^r < 2. The lower inset shows the same 
results normalized to the central NLO curve. The NNLO effects range from about 40% at pt ~ 50 
GeV to about 30% at px ~ 400 GeV. The NLO (NNLO) uncertainty ranges from about ±15% 
(±10%) at px ~ 50 GeV to about ±20% (±8%) at px ~ 400 GeV. We note that the NLO and 
NNLO bands do not overlap in the region where px < 300 GeV. This implies that, in this region 
of transverse momenta, the size of the band obtained through scale variations at NLO definitely 
underestimates the theoretical uncertainty. 

We now move on to the resummed results. In Fig.[4](a) the NLL±NLO spectrum is compared 
to the fixed-order NLO result and to the finite component of the resummed cross section (see 
Eq. (|4j) ) in the region between 0 and 80 GeV. As expected, the NLO diverges to ±oo as px —> 0, 
while the resummation provides a physically well behaved spectrum down to low values of px, 
which exhibits a kinematical peak at px ~ 4 GeV. The finite component contributes less than 1% 
in the peak region, where the result is dominated by resummation, and it increases to ~ 18% at 
px = 50 GeV. The lower inset shows the NLL±NLO result normalized to NLO. In Fig.[4](b) the 
region between 80 and 400 GeV is displayed. We see that even at large values of px the NLL±NLO 
resummed result does not match very well the fixed-order NLO result, with a difference of about 
5%. 


The analogous results at NNLL±NNLO are shown in Fig. [5} The NNLO has an unphysical 
(divergent) behaviour as px —>• 0, whereas the resummed spectrum is well behaved, with a slightly 
harder peak with respect to the NLL±NLO. The finite component contributes less than 1% in 
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(a) (b) 

Figure 4: The transverse-momentum spectrum of the W + W~ pair at NLL+NLO (a) 
in the low-p^ region and (b) at high transverse momenta. The NLL+NLO result (red, 
dashed) is compared to the fixed-order NLO prediction (grey, dash-dotted) and to the 
finite component of Eq. Q (magenta, dash-double dotted). The lower insets show the 
NLL+NLO result normalized to NLO. 




(a) (b) 

Figure 5: The transverse-momentum spectrum of the W + W~ pair at NNLL+NNLO (a) 
in the low -p? region and (b) at high transverse momenta. The NNLL+NNLO result (red, 
dashed) is compared to the fixed-order NNLO prediction (grey, dash-dotted) and to the 
finite component of Eq. Q (magenta, dash-double dotted). The lower insets show the 
NNLL+NNLO result normalized to NNLO. 
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p T (W + W') [fb] p T (W + W') [fb] 

(a) (b) 

Figure 6: W + W~ transverse-momentum distribution at the NLL+NLO (red, dashed) 
and NNLL+NNLO (blue, solid); thick lines: central scale choices; bands: uncertainty 
due to (a) /ip variation and (b) Pr variation; thin lines: borders of bands. 


the peak region, increasing to ~ 19% at p r = 50 GeV. Comparing the right panels of Fig. [4] and 
Fig.[5j we see that the quality of the matching at high p? is significantly improved when going 
from NLL+NLO to NNLL+NNLO, and we find that this behaviour is indeed preserved up to very 
high transverse momenta. The NNLL+NNLO result thus gives a prediction with uniform accuracy 
from small to very large transverse momenta and, in fact, provides a sufficiently large region where 
a hard switching to the fixed-order result is feasible. We point out that, thanks to our unitarity 
constraint, both at NLL+NLO and at NNLL+NNLO the integral of the resummed spectrum is in 
excellent agreement with the respective total cross sections; the differences are at the few-permille 
level. 


We now turn to the scale uncertainties of our resummed results. We start our discussion by 
separately considering factorization and renormalization scale variations. In Fig. [6] we compare 
the NLL+NLO (red, dashed) and NNLL+NNLO (blue, solid) predictions with their uncertainty 
bands from pp and Pr variations (left and right panel, respectively). In both cases, the bands 
are obtained by varying the factorization (renormalization) scale by a factor of two around its 
central value, while keeping the other scales at their default values. First of all, we notice that 
when going from NLL+NLO to NNLL+NNLO the pp spectrum becomes harder. Comparing with 
the results of Ref. 55 , where the NNLL resummation was implemented without C?(a|) matching, 


we see that the increased hardness of the pr spectrum is a combined effect of both features, i.e. 
NNLL resummation and NNLO matching at high pp. 


We note that neither in the case of the factorization scale, nor in the case of the renormalization 
scale, the NLL+NLO and NNLL+NNLO bands overlap. Actually, in the case of the factorization 
scale, there is no reduction in scale dependence when going from NLL+NLO to NNLL+NNLO, 
and the uncertainty slightly increases with the perturbative order, even if it is always well below 
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p T (W + W') [fb] p T (W + W') [fb] 

(a) (b) 

Figure 7: W + W~ transverse-momentum distribution at the NLL+NLO (red, dashed) 
and NNLL+NNLO (blue, solid); thick lines: central scale choices; bands: uncertainty 
due to (left) Pf,Pr variation and (right) Q variation; thin lines: borders of bands. 


10%, except at very low pp. The renormalization scale dependence instead exhibits the expected 
reduction when going from NLL+NLO to NNLL+NNLO. 

In Fig. [7] we present our resummed predictions with uncertainty bands obtained from simultane¬ 
ous variations of pf and Pr (left panel) and the variation of the resummation scale Q (right panel). 
In the left panel the uncertainty bands are obtained by varying pp and pn as in Fig.[3j In the right 
panel the resummation scale is varied in the range m^/2 < Q < 2 rriw As in Fig.[6]we see that 
the uncertainty bands do not overlap. The uncertainty from pp and Pr variations is ±10 — 15% 
at NLL+NLO and is reduced to 8 — 10% at NNLL+NNLO. At NLL+NLO the resummation 
scale uncertainty is generally about ±15% except in the region of pp ~ 10 GeV, where it shrinks 
to smaller values. We find that at NNLL+NNLO the resummation scale uncertainty is reduced 
roughly by a factor of two in the region of transverse momenta considered in the figure. 

In Figs.[8]and[9]we show our reference resummed prediction for W + W~ and Z Z. respectively, 
with an estimate of their full perturbative uncertainty. I 11 order to obtain a combined uncertainty 
from pp-, Pr and Q variations, we follow Ref. 1901 and independently vary pp, Pr and Q in the 
ranges my < {pfjPr} < 4 my and my/2 < Q < 2 my with the constraints 0.5 < Pf/pr < 2 and 
0.5 < Q/pr < 2. We recall that the constraint on Pf/pr, which is the same as applied in Fig. [3] 
and Fig. [7] (left), has the purpose of avoiding large logarithmic contributions from the evolution of 
parton densities. Analogously, the constraint on Q/ Pr avoids large logarithmic contributions in 
the expansion of the Sudakov form factor. 

For W + W~ production the perturbative uncertainty at NNLL+NNLO (NLL+NLO) is about 
±8% (±12%) at the peak, it decreases to about ±3% (±5%) at pt = 20 GeV, and it increases 
again to ±10% (±15%) at pp = 200 GeV. In the high-py region, the difference between the 
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p T (W + W') [fb] 



p T (W + W') [fb] 


(a) (b) 

Figure 8: (a) Transverse-momentum distribution of the W + W~ pair at NLL-f-NLO 

(red, dashed) and NNLL+NNLO (blue, solid); thick lines: central scale choices; bands: 
uncertainty from fj,p, fiR and Q variations obtained as described in the text; thin lines: 
borders of bands, (b) detail of the low -px region. 



Pt(ZZ) [fb] 



Pt(ZZ) [fb] 


(b) 


Figure 9: 


Same as Fig.[8j but for ZZ. 
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NNLL+NNLO and NLL+NLO predictions is driven by the NNLO effects, which increase the NLO 
result by about 30%. 

For ZZ production the uncertainties have essentially the same pattern in the small- and 
intermediate-pr region, while at high px they are larger than for W + W~ production, reaching 
about ±17% at NNLL±NNLO for px = 200 GeV. We have checked that this effect is entirely 
driven by the resummation-scale dependence. As previously pointed out, this behaviour is not 
particularly worrying since, in the large-pr region, the resummed results should be replaced by 
the corresponding fixed-order prediction. Also in the ZZ case the large enhancement of the 
NNLL±NNLO distribution in the liigh-p T tail stems from the fixed-order cross section. 


3.3 Rapidity dependence of the transverse-momentum distribution 



(b) 



Figure 10: (a) Shapes of the W + W~ transverse-momentum distribution differential in the 
rapidity of the W + W~ pair at the NNLL±NNLO for \y\ < 0.5 (red, solid), 0.5 < \y\ < 1 
(blue, dashed), 1 < \y\ < 2 (black, dotted), 2 < \y\ < 3 (magenta, dash-dotted), 3 < \y\ 
(orange, double-dash dotted); and (b) the shape-ratio with respect to the inclusive result. 


So far, we only considered px spectra for on-shell W + W~ and ZZ production that are inclusive 
in the kinematics of the vector-boson pair. Our numerical program, however, allows us to compute 
arbitrary observables that are differential with respect to the VV' phase space. 

In the following we study the behaviour of the transverse-momentum spectrum in different 
rapidity regions of the vector-boson pair. In Fig.[l0]we study the shape of the NNLL±NNLO 
transverse-momentum distribution, i.e. normalized such that its integral yields one, for \y\ < 0.5 
(red, solid), 0.5 < \y\ < 1 (blue, dashed), 1 < \y\ < 2 (black, dotted), 2 < \y\ < 3 (magenta, 
dash-dotted) and 3 < \y\ (orange, dash-double dotted). The right panel shows the same results 
normalized to the fully inclusive distribution. We clearly see that the px shapes become softer as 
the rapidity increases. In the central region (\y\ < 2) the distributions are still quite insensitive 
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to the specific value of the rapidity and only slightly harder than the inclusive spectrum. In the 
forward rapidity region, on the other hand, the shapes become increasingly softer. 


The observed pattern can be understood in the following way: rapidity and transverse momen¬ 
tum are two not completely independent phase-space variables. Indeed, they affect their mutual 
upper integration bounds. At higher rapidities the kinematically allowed range of transverse 
momenta is reduced: this squeezes the px spectrum which consequently becomes softer. This effect 
has been observed also in previous studies in the case of Higgs boson production 
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3.4 The W + W cross section and pp-ve to efficiencies 



NLL+NLO 

NNLL+NNLO 

NLO 

NNLO 

approx. NNLL+NLO 



Figure If: Veto efficiency for the transverse momentum of the W + W~ pair at various 
orders: NLL+NLO (red, dashed), NNLL+NNLO (blue, solid), NLO (grey, dash-dotted), 
NNLO (black, dotted), approximate NNLL+NLO (magenta, dash-double dotted); thick 
lines: central scale choices; bands: uncertainty due to combined scale variations; thin 
lines: borders of bands. 


The excess in the W + W~ production cross section measured by ATLAS 51 with respect to 
the SM prediction has drawn a lot of attention to the W + W~ 
appears in many new physics scenarios [53] 


process, since the W + W~ signature 
The inclusion of the recently computed NNLO 
considerably reduces the significance of the excess. However, particular attention 
must be payed to the modelling of the jet veto 52,58,93 when extrapolating from the fiducial 


corrections 36 
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region to obtain the inclusive cross section. Effects of jet veto resummation have been considered 
in Refs. 


94,95 , though still matching to the fixed-order 0(a$) result. 


In this paper we are dealing with transverse-momentum spectra, and we perform a resummation 
on a different variable with respect to the jet p t- However, the vector-boson pair px and the 
jet p T are clearly related variables (actually, at O(ag) they indeed coincide). We will therefore 
study the p^-veto efficiency in W + W~ production at different orders in resummed and fixed-order 
perturbation theory. We define the p^-veto efficiency as 


e(Px °) = a (PT < Pt °)/otot • 


(13) 


In Fig. 11 we show e(p™ to ) at the NNLL+NNLO (blue, solid), approximate NNLL+NLO (magenta, 


dash-double dotted), NLL+NLO (red, dashed), NNLO (black, dotted) and NLO (grey, dash-dotted). 
The lower inset shows the same curves normalized to our reference prediction at NNLL+NNLO. Our 
approximate NNLL+NLO is obtained by simply adding the g < A function in the Sudakov exponent 


in Eq. (10) at NLL+NLO, and corresponds to the approximation considered in Refs. 55,58 


For reference, the corresponding numerical values of the efficiencies are given in Tab. [I] for 
Pt = 5-40 GeV. The uncertainty bands are obtained by a combined variation of resummation, 
factorization and renormalization scales as in Fig. [8j The first thing we observe is that the NLO 
result appears to be well above the others and cannot be really considered a reliable prediction for 
the efficiency. This is because it is essentially a LO prediction at finite values of p™ to . We also note 
that in the small-pr region (say below pt ~ 10 GeV) the fixed-order NLO and NNLO predictions 
diverge and cannot be trusted. Comparing further the fixed-order results among each other and 
the resummed results among each other, we observe that higher-order corrections in fixed-order 
and resummed perturbation theory reduce the p^-veto efficiency. 


Both effects can be easily understood in the light of the results presented up to now. As seen 
in Fig.[3| the inclusion of the NNLO corrections make the Pt distribution harder. Furthermore, 
resummation effects generally harden the spectrum. A qualitatively similar result is obtained when 
going from NLL+NLO to NNLL+NNLO (see Fig.[8]). 


It is interesting to compare the approximated NNLL+NLO result with the NNLO and 
NNLL+NNLO predictions. For values of p™ to ~ 25 — 30 GeV we see that the approximated 
result is in between the NNLO one and our best NNLL+NNLO prediction. This means that 
the effect of NNLL resummation obtained by the inclusion of the r/ 3) function in the Sudakov 


exponent in Eq. (10) is quantitatively important. Nonetheless, the efficiency obtained within this 


approximation is still about 5% higher than the NNLL+NNLO prediction. We also notice that in 
this region of p™ to , the NNLO and NLL+NLO results differ by less than 1%. 

Comparing the NNLL+NNLO and NLL+NLO results, we find that they are compatible within 
the corresponding uncertainties. 


We add few comments on the recent measurement of the W + W 
CMS collaboration 


52 


cross section carried out by the 
The result shows good agreement with the NNLO prediction of Ref. 36 


The corresponding analysis, however, is based on a reweighting procedure of the px spectrum of 
the W + W~pair. The events generated with POWHEG 196 plus Pythia6 |97| were reweighted 
by using the calculation of Ref. j58], which corresponds to our NNLL+NLO approximation, and 
includes neither the second-order hard-collinear coefficient in Eq. (fob, nor the NNLO 


matching. The results in Fig. 0 show that the NNLL+NNLO p^-veto efficiency is lower than 
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77.2 
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74.0 

80.2 
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70 7+2.0% 

76.4 

77 q+1.3% 

1 / -°-3.6% 

77.1 

82.6 


Table 1: Predictions for the p^-veto efficiency (in percent) at various perturbative orders. 


the efficiency obtained with the approximated NNLL+NLO calculation. As a consequence, a 
reweighting to the full NNLL+NNLO prediction for the W + W~ spectrum would most likely lead 
to a decrease of the jet-veto efficiency. 


4 Summary 


In this paper we have studied the transverse-momentum distribution of vector-boson pairs in 
hadronic collisions. We presented a computation of the px spectrum in which the logarithmically 
enhanced contributions at small pt are resummed up to NNLL accuracy and the ensuing result 
is combined with state-of-the-art 0(ag) (NNLO) predictions valid at large px- We presented 
numerical results for W + W~ and ZZ production at the LHC together with a study of their 
perturbative uncertainties. 

We found that, up to relatively large transverse momenta, when scale variations are studied 
around the fixed resummation scale Q = my we obtain results which are fully consistent with those 
obtained using the dynamical choice Q = Myy /2. At large transverse momenta the fixed-order 
result in the tail of the distribution is nicely recovered with the fixed-scale choice. Our new 
NNLL+NNLO results significantly reduce the theoretical uncertainties obtained through scale 
variations compared to lower orders in both the peak region and the tail of the distribution. 

We have also studied the rapidity dependence of the resummed transverse-momentum distribu¬ 
tion. The rapidity dependence at NNLL+NNLO is quite flat in the central region (|p| < 2), but 
signals a substantially softer spectrum in the forward region. Due to phase-space suppression, the 
effect on the inclusive transverse-momentum distribution is very moderate though. 

Finally, we have studied the p^-veto efficiency at different orders in resummed and hxed- 


18 



order perturbation theory. Both NNLL resummation and the NNLO effects turned out to be 
important to obtain an accurate prediction for this quantity. We observed that the veto efficiency 
at NNLL+NNLO is ~ 5% lower (~ 3% in absolute terms) with respect to the approximate 
NNLL+NLO calculation used in the CMS analysis of Ref. 1521. This result suggests that our 
NNLL+NNLO predictions will be useful to validate the transverse-momentum spectra obtained 
from Monte Carlo event generators, similarly to what was done for the NNLL+NNLO calculation 
of Ref. 


59 in the case of Higgs boson production. 


In this paper we considered the px spectrum of stable vector-boson pairs, 
loop helicity amplitudes for qq —> VV' —> 4 leptons 137,38 


Exploiting the two- 
will allow us to extend the calculation 
to include the leptonic decay of the vector bosons and off-shell effects. The computation of the 
transverse-momentum spectrum with realistic experimental cuts will then become possible. 
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